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Abstract — The parameters affecting particle injection in a high pressure powder fuel injection system were verified by 
using computational particle fluid dynamics (CPFD). The particles were coal (d p =175 pm, p s = 1350 kg/m 3 ), and the 
simulation was performed by changing the clearance between the roller and the hose and by changing the total pressure 
drop. Also, the wall erosion of the hose was confirmed by changing the clearance. As the clearance decreases from 10 mm to 
1 mm when the total pressure difference is 10 bar, it is confirmed that the injection of the particles is increased two-fold, but 
the power of the motor should be increased to 1.45 times. When the clearance is 5mm, as the total pressure drop is increased 
from 10 bar to 30 bar, the particle injection rate decreases 0.33-fold, but the motor power must be increased 3-fold. Also, it 
was confirmed that the wall erosion of the hose was large when the clearance was 7 mm or more. 
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I. Introduction 

Many devices in the chemical process operate at high pressure rather than atmospheric pressure. In the case of the gas phase 
reaction, a high pressure induces large gas density. Therefore, a greater amount of gas can be added, which leads to an 
increase in the production amount. Recently, the Integrated Gasification Combined Cycle (IGCC), which has been attracting 
attention due to its low carbon policy, has shown high efficiency as well as low SOx, NOx, and C02 emissions compared to 
coal-fired power generation. However, the disadvantage of IGCC is its high cost compared to existing coal-fired power 
plants. Coal injection contributes about 43%-45% of the high cost of IGCC [1]. In the IGCC, it is necessary to inject coal 
from atmospheric pressure to high pressure; however, a general injection device cannot inject particles from atmospheric 
pressure to high pressure. This is due to the difficulty of injecting particles from atmospheric pressure to high pressure 
because the gas flows from a high pressure to a low pressure. In order to solve the problem of the backflow of gas, a lock 
hopper method is generally used when injecting particles from atmospheric pressure to high pressure. The advantage of the 
lock hopper method is that it can overcome a large pressure difference and can transfer the particles. However, the lock 
hopper method has disadvantages, since it is difficult to continuously inject particles because the pressurization and 
decompression must be repeatedly performed. Also, the lock hopper method requires equipment that is considerably large 
and relatively expensive in terms of construction and operating costs [2-5]. For example, in Shell’s IGCC plant in 
Buggenum, which the existing lock hopper method is used, the cost of the coal injection part is estimated to be about 0.1 
billion dollar [1]. Accordingly, a high-pressure powder fuel injection device has been developed to replace the lock hopper 
method injection device [1]. Stamet and Rocketdyne are developing the posimetric pump and the dry solid pump, 
respectively [1]. 

In this study, the experiment was carried out for the case of an atmospheric pressure condition, and the total pressure drop is 
10 bar. However, it is difficult to find the variables that may affect the higher pressures and injection rate. In addition, to 
check the wall erosion of the hose, the hose must be operated for a long time; therefore, the wall erosion needs to be 
predicted for validity analysis. Therefore, in this study, we predicted the solid transfer characteristics of a high-pressure 
powder injector by using computational particle fluid dynamics (CPFD) method using the multi-phase particle-in-cell (MP- 
PIC) method. 


II. Experiment and simulation 

2.1 Experimental 

Experiments were conducted using coal for verification of the simulation data. The specifications of the experimental 
apparatus are shown in Figure 1. A compression roller is installed in the chain, and the rotational force of the electric motor 
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is transmitted to the chain to rotate the compression roller. The upper hopper is filled with coal and the lower roller is 
compressed. At the same time, the high-pressure hose is pressed to transfer the coal in the high pressure hose to the lower 
hopper by the pushing force of the compression roller. The clearance between the hose and the roller was 10 mm. The speed 
of the roller was changed to0-0.12m/s using a motor inverter and the pressure of the lower hopper was changed to 10 bar 
at atmospheric pressure. 

The injected rate was calculated using the mass value of the sample measured, using a load cell for 10 minutes under each 
condition. 




Figure I.Schematic diagram of 

EXPERIMENT 



Figure 2.3D dimension of hoses that are pushed 

BY ROLLER. (A): REALIZATION MODEL, (B): 
SIMPLIFICATION MODEL 
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2.2 Governing equation 

The CPFD method calculates the momentum of particles and fluids in three dimensions. The fluid is calculated using the 
Navier-Stokes equation and the particle is calculated using the MP-PIC numerical method. The phases of the fluid and the 
solid are considered together by calculating the drag force between the phases. The phase of the fluid is solved by the 
governing equations. 

£M/) + V-(p/V/)= s / (1) 

For incompressible fluids, the momentum equation is expressed as (2). 

y t ( PfOfU f ) + V ■ (/ o f O f v f v f ) = -Vp - F + p f O f g + V6 f x f (2) 

where p f is the fluid density, 0 f is the volume fraction, v f is the fluid velocity, and S f is related to the fluid mass. x f is the fluid 
stress tensor, g is the gravity acceleration, and F is the ratio of momentum exchange in solid and fluid. F is calculated using 
equation (3). 

F = Iff fVp Pp [d p (v f - v p )- 1 - Vp] dV p dp p dv p (3) 


where v p is the particle velocity, p p is the particle density, D p is the drag function depending on the particle position, and f is a 
probability density function calculated using the Liouville equation. The particle acceleration is calculated using the 
Largrangian method. 


dvp 

dt 


= D p (v f -v p ) 


—VP + g-^—VT v 

pp ** Q v p v v 


(4) 


In equation (4), the value of D p is calculated using the drag model. The Wen-Yu Ergun model of Gidaspow is used, which 
combines the Wen and Yu model, known to be well suited for the dilute system. Also, the Ergun model is used, which is 
known to fit well in the packing system. Equation (5) represents the expression of the drag model. 

( D j, Op 0.75 O^p 

- °i) UCZZJ + • 0 ' 75f '» £ £ (5) 

D 2 Op 0.85 O^p 


where 0 p is the particle volume fraction and 0 cp is the volume fraction when the particles are packed at maximum. is 
calculated using the Wen and Yu drag model. 



PfWf-u p \ 
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( 6 ) 


where C d is the drag coefficient and is calculated using equation (7). 

Re < 0.5 
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The Ergun equation is expressed as equation (8). 
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0.5 < Re < 1000 
Re > 1000 


(7) 


( 8 ) 


2.3 Modeling of apparatus 

A schematic diagram of the 3D dimension for simulating the experimental setup used in this study with CPFD is shown in 
Figure2. The diameter of the hose is 0.04 m, the height is 1.45 m, and the thickness is 0.0125 m. The diameter of the roller 
pressing the hose is 0.16 m, and the distances between the first roller and the next roller are arranged at 0.3 m intervals. 
Based on the above information, a realization model is used when the roller presses the hose, and a simplification model is 
used that simplifies the cylindrical hose to a rectangular parallelepiped. In the case of the simplification model, the horizontal 
and vertical lengths of the rectangular parallelepiped are adjusted. As a result, the parallelepiped was simplified to a hose 
with a width of 0.035 m, a length of 0.035 m, and a height of 1.45 m. 


Page | 50 




International Journal of Engineering Research & Science (IJOER) 


ISSN: [2395-6992] 


[Vol-4, Issue-5, May- 2018] 


2.4 Simulation condition 

The average particle size of the particles used in the simulation was 175 pm, the density was 1350 kg/m 3 , and air was used as 
the gas. The volume fraction of the particles was set at 0.6. The variables considered to be significant in the simulation are 
the clearance and total pressure difference between the inlet and outlet. The clearance was changed to 1-10 mm and the 
pressure drop was changed to 10-30 bar. 

III. Results and Discussion 

3.1 Experiment result 

When the clearance between the hose and the roller is 10 mm using coal, the injection rate according to the roller speed is 
confirmed when the total pressure difference is 0 bar and 10 bar. As shown in Figure 3, the injection rate increases linearly 
with the roller speed at each pressure condition. However, it was confirmed that the injection rate decreased more when the 
total pressure drop was 10 bar than when the pressure drop was 0 bar. This is considered to be because, as the pressure at the 
lower end increases, the gas at the lower end is forced to flow backward due to the pressure difference, and the particles must 
descend to overcome the force. 



Figure 3. Roller speed versus solid feed rate 



Pressure, [MPa] 

Figure 4. Pressure profile of realization and 

SIMPLIFICATION 


3.2 Comparison of realization and simplified model 

In order to verify the validity of the simplified model, a simulation of the realization model is shown in Figure 4. By 
comparing the pressure drops among the rollers, an error of about 7% was found in the first roller at the bottom, 3% in the 
second roller, and 13% in the third roller. Although the contact areas of the rollers are the same, because the cylindrical roller 
is simplified into a rectangular parallelepiped, the difference is caused by the difference in height among the rollers. In 
addition, unlike the simplified model, the rubber hose is not pressed radially; rather, the pressure drop among the rollers 
differs from the actual model because the roller hose is gradually depressed according to the height of the roller. Likewise, 
the rubber hose is simplified to a rectangular parallelepiped in a simplified model, but the actual model is a cylinder. In the 
simplified model, the area pressed by the roller is constant, but in the actual model, the section which becomes extremely thin 
is generated. Therefore, in the simplified model, the pressure drop applied to the pressed part is constant according to the 
height, but in the actual model, the pressure drop according to the height increases in the narrower part. Increasing the 
pressure drop at the narrower part means that the force applied to the narrower part increases and eventually the entire the 
roller is not subjected to the same force but the part where the hose is pressed deeply under goes a larger force. However, the 
main purpose of this study is not to confirm the absolute value of the pressure drop, but to observe the change of the pressure 
drop according to the clearance or the tendency to change according to the total pressure. In addition, because the difference 
among the values is within 13% maximum, it can be confirmed that the simplified model can explain relation pressure drop 
and clearance. Consequently, the simulations were carried out using a simplified model. 

3.3 Pressure profile according to clearance 

Figure 5 shows the simulation results according to clearance and total pressure drop. Clearance refers to the distance 
between the roller and the hose at the roller contact distance. As can be seen in the graph, as the clearance is decreased, the 
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difference in pressure between the lower part and the upper part of the roller increases. On the other hand, the difference in 
pressure at the portion where the roller is not pressed decreased as the clearance decreased. When the clearance was 1mm, 
the pressure drop from the lower end of the hose to the lower end of the first roller at the lower end was 4.8 kPa and the 
pressure drop at the first lower roller was 207.7 kPa. When the clearance was 10 mm, the pressure drop from the lower end 
hose to the lower part of the first roller at the lower end was 71.4 kPa and the pressure drop at the first lower roller was 150.9 
kPa. As the pressure drop at the portion where the roller is not pressed is increased, the possibility that the high-pressure gas 
at the lower end flows backward increases. As the width of the clearance is reduced, the stable particle injection becomes 
more advantageous, but as the width of the clearance is reduced, the pressure drop on the roller increases. The pushing force 
of the roller is proportional to the pressure drop applied to the roller, which is 1.38 times greater than 1 mm when the 
clearance is 10 mm. In other words, when the clearance is 10 mm, the force applied to the roller is greater as the pressure 
drop is larger at 10 mm than the 1 mm clearance; therefore, a greater force is required to drive the roller. By comparing the 
pressure drops applied to roller when clearance is 1 mm, it was confirmed that a high pressure drop was applied when 
position of roller goes up, 207.7 KPa in the lower part, 272.7 KPa in the middle part, and 490.9 KPa in the upper part. That 
is, the upper roller receives more force because of the relationship between the volume and the pressure of the gas. The 
volume of the gas increases as the pressure decreases, and the velocity of the gas increases due to the increased volume. As 
the velocity of the gas increases, the energy of the gas increases and it is confirmed that the upper roller undergoes a greater 
pressure drop. Comparing the pressure drop of the roller due to the pressure drop between the upper and lower parts, the sum 
of the pressure drop across the three rollers is 9.7 bar when the total pressure drop is 10 bar and clearance is 1 mm, and is 
29.3 bar when the total pressure drop is 30 bar and the clearance is 1 mm. The roller pressure was3.02 times. For a total 
pressure drop of 10 bar and a clearance of 10 mm, a pressure differential of 3.08 times greater was loaded to the roller. As the 
clearance increases, the pressure drop increase of the roller decreases as the total pressure increases. On the other hand, when 
the clearance is 1 mm, the pressure drop increases 2.43 times, while it increases 2.83 times at 10 mm. As the total pressure 
drop increases, the possibility of backflow of the gas increases, but if the clearance is narrow, the rate of increase of the 
pressure drop at the portion where the roller is not pressed is reduced. Therefore, when the clearance is narrow, it is relatively 
easy to inject particles at high pressure. 



Pressure, [Pa] 

Figure 5. Pressure profile with varied 

CLEARANCE AND TOTAL PRESSURE DROP 



Clearance, [mm] 

Figure 6. Reverse gas flow rate with 

CHANGING CLEARANCE AND TOTAL PRESSURE 


3.4 Backflow of gas 

Figure 6shows the amount of backflow of gas at the top of the high-pressure powder fuel injector according to the clearance 
and total pressure. 

The gas backflow rate increased with increasing clearance and total pressure drop. The data of the gas backflow amount is 
plotted with the differential pressure data from the hose where the roller is not pressed. 

Outlet gas mass flow rate = 0.0166 [ ~~ ] * Pressure drop(MPa) - 0.001 (9) 
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That is, it can be seen that as the pressure drop of the portion of the hose where the roller is not pressed increases, the gas 
backflow increases. The differential pressure of the hose part where the roller is not pressed can be calculated by equation 
(9). 

Pressure drop[MPa] = Total pressure[MPa] * (0.0322 * Clearance[mm] — 0.0188) + 0.1224 (10) 

It can be seen that the pressure drop increases when the clearance increases or when the total pressure drop increases. As the 
back flow of the gas increases, the mass flux of the downward flow of the solid particles must be increased in order to 
prevent the back flow of the gas under the actual operating conditions. Determining the solid mass flux is related to the speed 
of the roller. As the speed of the roller increases, the amount of dragging and dropping of particles during the same time 
increases, and the solid mass flux value increases; this can prevent gas backflow. 

3.5 Wall erosion of hose 

The above sections describe the power and speed of the motor according to clearance and the total pressure. This section 
discusses the wall erosion of the hose. If a problem arises when the erosion of the hose becomes serious, shutdown then must 
be frequently performed while operating the process, which results in increasing costs. Therefore, it is very important to 
know the conditions that cause the erosion. Equation (11) is used to calculate the wall erosion and is based on the default 
value provided by CPFD. 

Wall erosion = f(0) * m s i s * v p 3 - s (11) 

Where m is the mass of the particle, v is the velocity of the particle, and f(0) is the value dependent on the angle. The wall 
erosion value is arelative value used to analyze the location of the wall erosion and the conditions under which further wall 
erosion would occur. Figure 7 shows the wall erosion profile versus height, changing the clearance. In all cases, the wall 
erosion of the hose occurred at the part where the roller was pressed. At the point where the roller begins to be pressed, the 
velocity of the gas increases, and the velocity of the particle increases, causing considerable wall erosion. Also, at the end of 
the roller, the accelerated particles dispersed and caused erosion on the hose. On the other hand, the wall erosion value was 
very large when the clearance was 7 mm or more. As can be seen in Figure 8, when the clearance is greater than 7 mm, the 
movement of the particles in the x-direction is freer and the particle x velocity value is increased compared to when the 
clearance is narrow. The increase in particle x velocity resulted in more erosion effects than with particle z velocity; this 
caused a greater amount of erosion with a clearance greater than 7 mm. 



Wall erosion 

Figure 7.Wall erosion versus axial height with changing clearance 
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Figure S.Particle velocity versus axial height 


3.6 Prediction of high pressure powder fuel injection system using CPFD 

Simulation shows the effect of clearance and total pressure difference. In the experiment mentioned in Section 3-1, the 
injection amount ( m s ) of the particles according to the speed (0024-0.12 m/s) of the roller at the atmospheric pressure is 
expressed by equation (13) at equation (12) and 10 bar. 

Solid mass flow rate((ra s ) [yj = 2049.3 x Roller speed + 2.5168 (12) 

Solid mass flow rate((m s l0 ^ ar ) [^] = 789.15 [^] x Roller speed [“] + 3.5896 (13) 

The slopes of equation (13), i.e., the weights of the particles per unit length, were compared using data obtained from CPFD 
to experiment. Assuming that the gas backflow is zero at normal pressure, the weight of particles per unit length is 2049.3 
kg/m and the weight of particles per unit weight is 789.15 kg/m when the backflow is 0.4256 kg/s. Based on this data, Table 
1 shows the predicted motor force based on the predicted value of clearance and the particle weight per unit weight for the 
total pressure drop and the pressure drop on the roller, with 10 bar and clearance 10 mm as references. If the roller is 
operated at 20 bar at intervals of 10 mm, then it is not operating normally because it has a value of -0.108.6 kg/m. It can be 
seen that it is possible to reduce the gap to 5 mm or less in order to obtain a similar injection amount. However, the decrease 
in clearance increases the pressure drop of the roller, which means that the power of the motor increases. For example, if the 
clearance is 5 mm at 20 bar, the motor power needs to increase 2.53 times more than the reference. 


Table 1 

Predicted value of mass per unit length and power ratio 


Pressure, [MPa] 

Clearance, [mm] 

Mass per unit length, [kg/m] 

power ratio, [-] 

0 

10 

2049.3 


1 

10 

789.1 

1.00 

1 

7 

1075.2 

1.14 

1 

5 

1265.8 

1.25 

1 

3 

1456.5 

1.35 

1 

1 

1647.2 

1.45 

2 

10 

-108.6 

2.05 

2 

7 

463.4 

2.32 

2 

5 

844.8 

2.53 

2 

3 

1226.2 

2.72 

2 

1 

1607.5 

2.91 

3 

10 

-1006.3 

3.08 

3 

7 

-148.3 

3.53 

3 

5 

423.8 

3.82 

3 

3 

995.8 

4.10 

3 

1 

1567.9 

4.37 
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IV. Conclusion 

CPFD simulations were carried out to identify the variables affecting the operation of the high-pressure powder fuel injection 
system and to understand the effect of the variables. It was found that as the clearance increased, the backflow of the gas also 
increased. However, as the clearance increased, the pressure drop on the roller decreased, so the load on the motor was 
reduced. On the other hand, CPFD simulation predicted that a clearance of 5mm or less should be maintained for the wall 
erosion of the hose, since a clearance of 7mm or more is expected to cause significant erosion. 
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Nomenclature 


dp 

Sauter mean diameter [□ m] 

m s 

Mass of solid [kg] 

Ps 

Solid density [kg/m 3 ] 

e 

Angle of collision [°] 

Pp 

Particle density [kg/m 3 ] 

s f 

Interpolation operator [-] 

pf 

Fluid density [kg/m 3 ] 

Tf 

Stress tensor of fluid [Pa] 

Qp 

Particle volume fraction [-] 

g 

Gravitational acceleration [m/s 2 ] 

0cp 

Close pack particle volume fraction [-] 

F 

Momentum exchange rate per volume between the 
fluid and solid phases [Pa*s/m] 

Of 

Fluid volume fraction [-] 

D P 

Interphase drag coefficient [m/s 2 ] 

Vf 

Fluid velocity [m/s] 

f 

Particle distribution function [-] 

v p 

Particle velocity [m/s] 

c d 

Drag coefficient [-] 

rh s 

Solid mass flow rate [kg/s] 
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